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Abstract: We study Gravitational Waves (GWs) in the context of Massive Gravity, an 
extension to General Relativity (GR) where the fluctuations of the metric have a nonzero 
mass, and specifically investigate the effect of the tensor modes on the Cosmic Microwave 
Background (CMB) anisotropies. We first study the time evolution of the tensor modes in 
Massive Gravity and show that there is a graviton mass limit mi = lQ~ m g ~ 10 -29 cm -1 , so 
that for masses m < mi the tensor perturbations in Massive Gravity are indistinguishable 
from the corresponding ones in GR. Also, we show that short wavelength massive modes 
behave almost indistinguishably from their massless counterparts. Later on, we show that 
massive gravitons with masses within the range m = 10 - cm - m = 10 _26 cm _1 would 
leave a clear signature on the lower multipoles (1 < 30) in the CMB anisotropy power 
spectrum. Hence, our results indicate that CMB anisotropies measurements might be 
decisive to show whether the tensor modes are massive or not. 
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1. Introduction and Summary 

Over the last years we have been witnessing a great evolution in the observation of the 
universe. An accurate statistical analysis of the anisotropies of the CMB using the five- 
year WMAP data [|J has shown that the most favored cosmological model to fit the data 
is the flat ACDM model, which includes not only the very well known baryonic matter, 
but the mysterious dark matter (DM) and dark energy (DE) as well. The introduction 
of dark matter does not require a modification of GR (despite some alternative models 
do); in the pure general-relativistic case, all we have to do is to add further terms to the 
energy-momentum tensor on the right-hand side of the GR field equations. The DM can- 
didates must pass some important tests Q to ensure that the resulting model is physically 
consistent. However, despite the good theoretical candidates available in f2|, there is no 
observational evidence so far to support any of them as the actual components of DM. 

However, the inclusion of DE is not that simple. Ordinary matter, either baryonic or 
dark, cannot accelerate the universe as shown by observations Q, 0]. This discovery posed 
one of the deepest questions in modern cosmology: is the universe accelerating due to a 
repulsive gravity (caused by the quantum energy of the vacuum, for example, as in the 
ACDM model), or does General Relativity (GR) break down on cosmological scales 0? 



Many attempts have been made addressing these two possible cases, either by introducing 
new features into GR, or by modifying it (see || for a review of the models for DE). 

The two key points introduced above are just to illustrate the need of studying al- 
ternative theories of gravity in parallel with the improvement of models in GR itself. As 
alternative theories we mean modifications of GR, and the simplest possibility for modify- 
ing GR is the introduction of a mass for the graviton. It was performed for the first time 
by the pioneering work of M. Fierz and W. Pauli where they considered a linearized 
field theory of spin- two massive particles. The Lorentz invariance of the Fierz-Pauli (FP) 
lagrangian yields a spin-two massive state with five polarization modes (states with helic- 
ities ±2, ±1 and 0), differing from GR where one finds only a spin-two state with the two 
tensor polarization modes (helicities ±2). Such extra degrees of freedom yield an additional 
contribution of one vector and one real scalar massless particles with helicities ±1 and 0, 
respectively. The scalar particle couples to the trace of the stress energy-momentum ten- 
sor, causing a discontinuity in the propagator when one switches from the massive to the 
massless regime. This is the so-called van Dam-Veltman-Zakharov (vDVZ) discontinuity 
11, ||, whose net effect for a theory of a massive spin-two graviton is catastrophic: it 
would not even pass the solar-system tests for a theory of gravity (the prediction of the 
angle concerning the bending of the light by the Sun, for example). 

However, in a full theory of gravity, we must consider nonlinear effects; the FP theory is 
valid only in the linear approximation. Nonlinear effects eliminate the vDVZ discontinuity 
in the classical level [pi]] , |lT[| , so that classically we may reconcile the massive theory with 
the GR predictions. 

However, at the quantum level, the nonlinear interactions appear at the loop diagrams, 
so that the theory becomes strongly coupled above the energy scale A = (m 4 Mp;) 1/ ' 5 , where 
m is the graviton mass and Mpi is the Planck mass [fL^] , [13]. For masses m ~ Hq, where 
Hq is the present-day value of the Hubble parameter, the energy scale A is too small, well 
below the expected value, A = {mMpi) 1 ! 2 . In brane-world models (fl4|, IS], [16], []17|]) 
a similar problem occurs: either they have ghosts JI8J , [|Ts|], |2C)[ ] and pi] , or are strong 
coupled at low energies 0, Q, @ and |22|. 



A great step forward was taken in the works [23] and [24|. In reference ]23j the authors 
proposed a consistent modification of gravity in the infrared as an analog of the Higgs 
mechanism in GR. In this model, Lorentz invariance is spontaneously broken and the 
graviton, as a result, acquires a mass. In reference [24] the author introduces a Lorentz- 
violating massive gravity model in which the vDVZ discontinuity, ghosts and the low strong 
coupling scale are absent. In reference [25| the author studies the most general Lorentz- 
violating gravitational theory with massive gravitons, showing that there is a number of 
different regions in the mass parameter space of this theory in which it can be described by 
a consistent low-energy effective theory without instabilities and the vDVZ discontinuity. 

Therefore, the theory of Massive Gravity, as developed in |24|] and p5|], gives rise 
to physical propagating modes, and is free of the pathologies mentioned above. It is a 
potential candidate to provide the proper answers to the open questions in cosmology 
as mentioned earlier. There is a number of works studying cosmology in the context of 
Massive Gravity |26[| , [27|, [28]; in the present work, we aim at extending this discussion by 



analyzing the anisotropics of the CMB induced by tensor perturbations in Massive Gravity. 
We have analyzed tensor and vector perturbations in theories of gravitation with massive 



gravitons in a previous work [29], and in this paper we focus specifically on the signatures 
of the massive tensor modes. Henceforth we shall use the terms graviton and tensor modes 
interchangeably. 

We know from a number of sources [ 30 ] , [31|, fljffi ], [33], that primordial GW might 
leave a signature in the anisotropies and polarization spectrum of CMB, generated by 
the influence of such GWs on the photon redshifts. So, if the gravitons do have a mass, 
we expect that they will leave a different signature on the CMB anisotropy spectrum. 
Therefore, the main goals of this present work are twofold: develop solutions for massive 
GWs and analyze their signatures on the CMB polarization spectrum. 

To this end, the present paper is organized as follows: in section |2| we review the 
basics of Massive Gravity and its cosmological tensor perturbations. In section |3| we start 
reviewing the solutions for primordial GW in GR, and right after it we present our results 
in Massive Gravity; furthermore, we compare the results of both theories for different 
graviton masses, and discuss the differences between them. In section |] we review the 
basics of radiative transfer in the presence of weak gravitational fields, deriving the relevant 
Boltzmann equations. In section || we derive the expressions for CMB anisotropies and 
polarization, with the corresponding correlation functions. In section ^ we discuss the 
solutions to the Volterra integral equation which provides the functions that enables us to 
evaluate the coefficients to the harmonic expansion the modes introduced in section In 
section ^ we apply all these theoretical tools to Massive Gravity, obtaining the individual 
power spectra for different wavenumbers and graviton masses, and we compare these results 
with the predictions of GR. At the end of this paper we discuss the obtained results and 
make the corresponding conclusions. 

2. A quick overview of Massive Gravity 

As we have pointed out in the Introduction, the key ingredient to construct a physically- 
consistent theory of gravitation with massive gravitons lies on the spontaneous violation 
of the Lorentz symmetry. As in the Higgs analog in the Standard Model of electroweak 



interactions, we introduce, following |25|] and [26], a set of four scalar Goldstone fields 
</>°(x), 4> l {x), such that the action for Massive Gravity is written as 

S = f d 4 x^g- [-M 2 Pl R + A 4 F(X, V\ W**) + C matter ] , (2.1) 

where the first term on the right-hand side represents the usual Einstein-Hilbert action, and 
F is an arbitrary function of the metric components, their derivatives, and the Goldstone 
fields. The lagrangian for ordinary matter, C ma aer, is assumed to be minimally coupled to 
the metric. The simplest way to combine the derivatives of the Goldstone fields to enter 
the argument of F is given by the set of scalar quantities 



X = A" Wd^ , 



(2.2) 



yi = A-^datfPdptj?, 



V i V j 
X 



(2.3) 
(2.4) 



where A is the parameter which characterizes the cutoff scale of the theory. The second term 
on the right-hand side of ( [2.1D is invariant under the spatial reparametrization symmetry 
x l (t) — > x l (t) + and rotations. 



We now introduce the "vacuum" solutions for the model ( 2.1 



9a(3 = a 2 r] a p, (jp = A 2 t, 



AV, 



(2.5) 



which corresponds to the flat FRW space; in the "unitary gauge" described by ( |2.5| ) the 
action will depend solely on the metric components. Now, in order to study linear cosmo- 
logical perturbations around a flat Friedmann-Robertson- Walker (FRW) space, we spon- 
taneously break the Lorentz symmetry of the model by fixing the Goldstone fields to the 
vacuum (|2.5[), so that the only remaining perturbations are given by 



9a/3 = a rj aj3 + 8g Q/3 , 



(2.6) 



where rj a p = diag{+, — , — , — }, a(rj) is the scale factor, and 5g a p is a metric perturbation 
whose components are given by pj], 



2 



2a 2 <p, 5g 0i = a 2 (Si - diB), 



5gtj = a [-hij - diQj - djQi + 2(^5^ - didjE)] 



(2.7) 
(2.8) 



where cp, ij), B, E are scalar fields, Qi and Si are vector fields, and hij is a tensor field. The 
constraints satisfied by the vector and tensor fields are PH , 34], 



hij' 3 = 0, h 1 



0, Q i 
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Now, in the unitary gauge ( |2.5| ) we expand \J —g + 6g, X(g+5g), V l (g+5g), W lJ (g+5g) 
and F(g + 5g) in powers of the metric perturbation 5g, and substitute these results into 



the massive term in (2.1), so that the lagrangian for the second-order perturbations reads 



[ m o^oo + 2m i s 9oi ~ mi s 9ij + ™%8gu8gjj - 2ml8g o5gu] , (2.10) 



where mo, m%, m%, m% and are parameters related to the function F and its derivatives, 
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(2.11) 

(2.12) 
(2.13) 
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where W = -l/2,5 i:j W ij and 



Fx 



_ OF 
~ dX' 

F w 5ij 



Fxx 



dF 



8 2 F 

dx 2 ' 



Fvv§i 



8 2 F 



Fxw^i 



8 2 F 



d 2 F 



FwwAjSki + Fww2(5ik5ji + $u5jk) 



(2.14) 

(2.15) 
(2.16) 



dW^dW kl 

(see Appendix A in references [26] and [28| for details). The spatial indices in ( 2.10[ ) are 
summed over and, as argued in the reference [PH , the mass parameters mi are proportional 
to some scale denoted by m. 



The Einstein equations for the model (2.1), with the Goldstone fields in the unitary 
gauge (|2.5| ), and metric ( |2,6| ) read (for computational details, see appendix A of the refer- 
ences p§ and |H), 



m 1 



in! + n 2 



a 



(Pm + Pcf, 



(Pm + 



Pa), 

><t> +PA), 



(2.17) 
(2.18) 



doiapFxX 1 / 2 ) = 0, (2.19) 

where TC = a' /a, p m and p m stand for the density and pressure for the ordinary matter 
respectively, and 

p^ = A 4 XFx, p^ = A 4 WF w , (2.20) 



Pa 



—F, 

2 ' 



PA 



—F. 
2 



(2.21) 



The prime represents a derivative with respect to the conformal time rj. 

Once we have established the dynamical equations for the background, let us now turn 



our attention to the metric perturbations (2.8). The steps toward obtaining the dynamical 
equations for the massive metric perturbations are quite similar to those referred in pi} ], 
and they can be found in details in the Appendix A 3 in |28[ ]; here we simply quote the 
results. Since in this paper we are interested solely in the tensor perturbations represented 
by the element hij in (|2.§|), we write its dynamical equation as [28] 



V 2 h i:j + W.hL + a 2 mlhij = 0. 



(2.22) 



Since we deal only with the mass throughout this paper, we henceforth drop the sub- 
script 2 and write it simply as m. 

To end this section let us discuss an important aspect concerning the mass parameters 
of Massive Gravity. As we have pointed out in the Introduction, there are regions in the 
mass parameter space in which this theory is free of ghosts and instabilities; this means 
that the mass parameters mo, mi, m,2, and 7714 cannot be chosen arbitrarily, but they 



have to satisfy some constraints [24], [25]. Since in this paper we deal only with the mass 
parameter ttt-2, there is a number of choices on these parameters in which the model is 



physically healthy; therefore, any of these choices would produce a physically acceptable 
theory. We simply assume that the mass parameters in our work are within the region in 
which the pathologies are absent. 

Specific restrictions on the function F are discussed in [|26[|. In this reference, the 
authors demonstrate the existence of a wide class of functions F for which expanding 
cosmological solutions are compatible with constant graviton masses and allow for the 
effective field theory description. Therefore, we may simply restrict F in such a way the 
mass ??i2 is constant along the story of the universe, which we assume to hold throughout 
this paper. 



3. Primordial Gravitational Waves in GR and Massive Gravity 
3.1 Primordial Gravitational Waves in GR 

In GR, a primordial GW is described by the cosmological tensor perturbation whose dy- 



namical evolution is governed by the equation 34] 



h'-j - V 2 hij + 2Hti i:j = 0. (3.1) 

Before going into the Fourier space to solve this equation, it is convenient to introduce a 
new parametrization into this model J35|]. Let us write down the present-day scale factor 
a(rjo) as a quantity with dimension of length; then, setting Rh = c/Hq as the Hubble 
radius, we then define 0(770) = 2Rh- Now, since the GW wavenumber k is very small for 
primordial GWs (in the frequency range which could produce a signature on CMB), with 
wavelength comparable to the present-day Hubble radius Rh, we introduce a dimensionless 
time-independent vector n which has the same direction of k, and whose modulus is exactly 
the proportionality factor between the modulus of k and Rjj: 

n = 2R H k. (3.2) 

Since the early cosmological perturbations are of quantum-mechanical origin, we construct 
the tensor hij as a quantum-mechanical operator, whose Fourier expansion is given by 



ha (t?,x) = ^|pp [^nK^e-X + 4;(n)C(r/)e- in - x a r n t], (3.3) 



^ipi f°° d 3 i 

n 

r=l,2 

where r stands for the polarization mode of the GW, is the GW polarization tensor, 
and Ipi is the Planck length. The annihilation and creation operators a T n and a n satisfy 
the well known commutation relations 

a r n ,a s j]=5rs5 {3 \n-n'), (3.4) 



and, for the vacuum state |0), 

a r a \0) = 0. 



(3.5) 



Now, substituting the expression (|3,3| ) into (|3.1| ), and redefining the conformal time deriva- 
tive as d/drj = (a/c)d/dt, we get the equation governing the dynamics of GR tensor modes 
(dropping the polarization index for a while) 



h!' n + 2Hh' n + n 2 h n = 0. (3.6) 

Then, defining the quantity n n (rj) = a(rj)h n (rj) fl3q] , we obtain an equation for a paramet- 
rically disturbed oscillator 



n 2 



a" 



Mn = 0. (3.7) 



a 

To solve the simplified equation ( |3.7[ ) we have to specify the scale factor. Since we are 
mostly interested in the time of recombination, we use the scale factor for a flat universe 
filled with radiation and matter, whose expression is |3S 



( , 2 V J +' 



where z e<? is the redshift associated with the epoch of radiation-matter equality, whose 
value is z eq ~ 3 x 10 3 , and the corresponding conformal instant r] eq is given by 



Veq = (V2- i) V 2 + ffg „ 7 6 x 10 -3 (3 Q) 

1 ~r ^eg 



It is easy to see that ( |3.8j ) reduces to a scale factor for a radiation-dominated and 
matter-dominated universe 

a(V) = , V, V < r] eq ; 

a(ry) = 2i?^ (77 + r/ e(? ) 2 , r\>r\ eq (3.10) 

respectively, so that it comprises the whole period we are interested in. 

Substituting the scale factor (|3.8|) into ( |3.7l ), we can obtain exact analytical solutions 



for the functions Univ) IeH- Following [33], we normalize the GW amplitudes h n (rj) in 
terms of its value at r\ r = 10 -6 (in terms of redshift, z r ~3x 10 7 ); the resulting numerical 
solutions are displayed in the figure ||. 

3.2 Primordial Gravitational Waves in Massive Gravity 

Once we have reviewed the properties and evolution of GW amplitudes in GR, let us now 
analyze the same issues in Massive Gravity. First of all, we treat the tensor perturbations 
in the massive case quantum-mechanically, so that the Fourier expansion for the massive 
tensor field hij (rj, r) has an analog expression as in 



+ e\™ )r *(n)h^^ V )e-™-*dt )r \ (3.11) 




Figure 1: The time evolution of the normalized GW amplitudes h n (r))/h n (ri r ). Compare with 
Figure 1 of |1. 



where the superscript (m) stands for massive, and E n denotes the energy of the mode n. 
Now, plugging fl3.11| ) into ( |2.22| ) we get 



/4 m) " + + (n 2 + m 2 a 2 ) fcM = 0. 



(3.12) 



In ( 3.12 ) we have dropped the GW polarization indices r because we are going to treat only 
the tensor modes in this paper. Using the same strategy as in the GR case, we introduce 
a function ^^(rj) = a(r])h^ '(77), so that equation ( |3.12 ) becomes 



/4r } " + 



2 , 2 2 

n + m a 



a 
a 



M = o. 



(3.13) 



As we have discussed in the section above, the scale factor (|3.8|) represents very well 
the periods of the universe considered in this paper, so that it makes sense to employ it in 
Massive Gravity as well, since we may expect that the contribution of massive gravitons to 
the expansion of the universe is negligible in its early epochs; then, as a first approximation, 
we may neglect the contribution of the components p^, equation ( 2.20 ), and p\, equation 
Q , in ( ET7D . 



Now, using the above arguments and consequently the scale factor (3.8), we can solve 
numerically equation ( |3.13| ) for different wavenumbers n and masses m. We choose the 
graviton masses m using the following argument: in GR, only GW with frequencies v 
within the range 10~ 15 Hz to 10~ ls Hz may leave a signature on CMB polarization 



these frequencies correspond to wavenumbers k within the range 10 -25 cm -1 (n ~ 5 x 10 3 ) 
to 10~ 28 cm -1 (n ~ 10). For Massive Gravity, we use the same values for k, but now 
we vary the frequencies in order to obtain constant nonzero graviton masses through the 
dispersion relation 

lo 2 = k 2 + m 2 , (3.14) 



which comes straight from ( [2.22| ), where now u = 2nv. As a result, we find that if the 
values of the mass m lie within the range 10" 66 - W~ 62 g, the corresponding frequencies 
have values very close to the expected in GR. In particular, we've found that if the graviton 
mass is m = 10~ 66 <? ~ 10~ 29 cm _1 , the behavior of the GWs in Massive Gravity is exactly 
the same of GWs in GR. Therefore, if the graviton mass is equal or less than the graviton 
mass limit mi = 10 -66 g, the effects of Massive Gravity are indistinguishable from GR. 
It is important to mention that there has been a lot of efforts to constrain the masses 



of the tensor modes over the past few decades. For instance, Goldhaber and Nieto [38] 
have found a limit m < 2.0 x lCP 62 g analyzing the motion of galaxies in clusters. Later on, 
Talmadge et al. [|39| studied the variations of Kepler's third law when compared with the 
orbits of Earth and Mars, and found a limit m < 7.68 x 10~ 55 g. Recently, Finn and Sutton 
flio[| calculated the decay of the orbital period of the binary pulsars PSR B1913+16 (Hulse 
and Taylor pulsar) and PSR B1534+12 due to emission of massive gravitons, and found 
m<1.4x 10~ 52 g. Cooray and Seto [41| investigated the variation of the speed of gravity 



when compared to the speed of light due to a massive tensor mode, and determined an 
upper limit of ~ 10~ 56 <? for its mass, by using the measurements of a sample of close white 
dwarf binaries detectable with the Laser Interferometer Space Antenna (LISA), together 
with a optical light curve data related to binary eclipses from meter-class telescopes for the 
same sample. A recent and comprehensive review of the methods to determine the bounds 
for the masses of gravitons and photons can be found in p^] , which we refer to for further 
details. 

Since we are interested in investigating signatures of massive gravitons, we shall con- 
sider only graviton masses higher than the limit m = 10 _66 g; the numerical solutions to 
the massive tensor perturbation equations ( p. 13 ) are depicted in the Figure || below. For 



sake of comparison we depict the general-relativistic GW amplitudes in each graph as well. 



We have used the same normalization as [33], and the plots start at rj = rj r = 10 . The 



mass m = 2.843 x 10 28 cm 1 correspond to m = 10 e5 g, and so forth. 

Let us now analyze in detail the behavior of massive gravitons in the light of equation 
(|3.13[) . In the very early universe, before the time of equality radiation-matter, the value 



of a{rj) is very low, and then the m a on the left-hand side of (|3.13[) can be dropped; 



therefore, we recover the characteristic tensor mode equation of GR, ( |3.7D , and the behavior 
of massless and massive gravitons are the same. On superhorizon scales, n <C a" /a, the 
resulting equation for the tensor modes is 



a" 



^" ~ = 0, (3.15) 

whose solution is given by fi^ = f( n ) a , which means that the tensor amplitudes are 
"frozen", no matter the gravitons are massless or not. This particular behavior can be 
clearly seen from figures |2] - @, where the amplitudes are constant for all the modes consid- 
ered prior to decoupling. 

However, as the universe evolves, the tensor modes "fall" into the horizon, so that 
their amplitudes are no longer constant; on subhorizon scales, n 3> a" /a, we can neglect 
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Figure 2: The time evolution of the normal- 
ized GW amplitudes for n = 5. 
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Figure 4: For n = 50. 
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Figure 3: For n = 10. 
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Figure 5: For n = 100. 
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the effect of the term a" /a, so that we are left with 



M" + [ n 2 + m 2 a 2] ^(m) = (3 16) 



On subhorizon scales the massive term becomes dominant over low values of n, so that it 
"enforces" the tensor modes to fall into the horizon earlier than in the massless case. It is 



clear from equation ( 3.16 ) that the heavier the gravitons, the earlier their modes fall into 
the horizon. 

This effect can be clearly seen in the figures ^, || and |4], where the n values are suffi- 
ciently low to account for this effect. However, for larger values of n, this effect weakens, 
since n 2 3> m 2 a 2 in the time of decoupling, and the massive term will be predominant only 
for low redshifts, as can be seen in figures || |6| and ^. 

In particular, for low n (corresponding to tensor modes with long wavelengths), its 
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Figure 6: For n — 500. 



Figure 7: For n = 1000. 
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Figure 8: The "tail of figure [| zoomed in, showing the phase difference in the tensor modes at 
very low redshifts for both massless and massive gravitons. 



constant contribution to ( 3.16| ) can be completely neglected, so that we are left with 

M M = 0, (3.17) 



2 , 2 2 a 

n + m a 



and then the oscillatory behavior is strikingly different from the massless shown 
in figures ||, ||. For higher n (that is, tensor modes with short wavelengths), tough, this 
effect is not so strong, but induces a slight phase difference in the oscillatory behavior of 
the tensor modes. Such phase difference is stronger for higher masses; as an example of it, 
we have zoomed in the "tail" of figure || to show this fact. This is presented in figure ||. 

Hence, from this analysis we may conclude that the tensor modes of Massive Gravity 
behave similarly to the massless modes of GR, but the heavier the tensor modes are, the 
more distinct are their physical evolution if compared to the massless modes. Nevertheless, 
if massive gravitons do exist, they likely have left a signature on some physical observable; 



then, by comparing the predicted signatures of the massless and the massive modes with the 
observed ones, one should be able to determine whether they possess or not a nonzero mass. 
In our view, the best laboratory to test this assumption is the anisotropies measurement 
of CMB, for reasons that will become clear in section [?|. Before doing so, though, we take 
some time to review the basic aspects of the theory of CMB anisotropies and polarization. 



4. The Radiative Transfer Equation in the presence of Weak Gravitational 
Fields - an overview 

The first account of the effect of primordial GWs in polarizing the CMB photons was 
introduced in a seminal paper by Polnarev (3(J. To begin with, let us consider a given 
beam of radiation characterized by its Stokes parameters Q] {/, Q, U, V}, where I is the 
total intensity of the wave, the parameters Q and U measure the linear polarization of 
the wave, and V measures its circular polarization. They are integrated over all radiation 
frequencies, so that there is a set of Stokes parameters for each monochromatic component 
wave of the radiation beam {I(y,9,<p),Q(v, 9, tp), U(y, 9, tp),V(u, 9, tp)}. Associated with 
the Stokes parameters are the components of the photon distribution function n, which 



can be cast in a symbolic vector of the form 3C], 

2 fi + Q 



- - 1 c 

n ~ 2/^3 



I-Q 

-2U 



(4.1) 



Now, the transfer equation for the photon distribution functions subject to a weak GW-field 
is given by 

C[n], (4.2) 



dh „ j <9n <9n dv 
dr/ dx % dv air] 



where e l is the unit vector along the photon geodesic, and C [n] is the scattering term given 
by 



C[n] = -a T N e a(rj) \ n(Tj,T,i/,fj,,<p) 



dfi'dip P (pt, <p, fjf, <p') h(7], r, i/, fjf, ip') \ , (4.3) 



47T J_i 

P (n, tp, ip 1 ) is the scattering matrix given by [|^] 

p (», tp, J, tp') = q|p° [fi, fjf) + Vi-^Vi-^ 1 (m, ^ a*', + p 2 (m, ^ f') } . 

(4.4) 



where 



Q 



pu 



(\ o o\ 
10 
10 
\0 2/ 

f 2(1- fi 2 ){l-fi' 2 )+ fi 2 fi' 2 fi 2 \ 

1 



fifjf J 



(4.5) 



fi 



(4.6) 
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o 
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P 1 



(4.7) 



( 4/z// cos tj) — 2/isin^ \ 



2 // sin tj) cos?/' 

\ cos V J 

( /U 2 // 2 cos 2ip — ^ 2 cos2V> —/i 2 ^' sin 2V' 0\ 
— /x' 2 cos 2-0 cos 2^ fj,' sin 2^ 
/i/Lt' 2 sin 2ip —[A sin 2-0 pp! cos 2-0 
V o 0/ 

or is the Thomson scattering cross-section, N e (rj) is the number of free electrons in the 
unit comoving volume, fi = cos 6, and we have defined ip := ip — ip 1 [30]. 

In order to get the expression of the Boltzmann equation for our problem, let us 
decompose the vector n into its zeroth-order contribution, n(°) i. e., in the absence of GW, 
and its first-order correction n^, 



(4. 



n = n(°)+n( 1 ), 



(4.9) 



where is the blackbody radiation function 



n(°V) 



1 



(4.10) 



and To ~ 2.725 K is the present-day value of the CMB temperature. Since equation ( p~2; ) 
is linear, we can expand in the same way as we did in ( [3.3D , 



J-oo V2n r ^ x I 



(2^) 3 / 2 



m-x-rf 



(4.11) 



now, introducing the basis for Thomson scattering [30], 



OrC/i, ip) = ~ (1 - ^ 2 ) e ±2 ^ u, 6 r ( M; = i 



/ (i + r 



(4.12) 



where r = 1 corresponds to a left-hand polarization, and r = 2 to the right-hand one, we 
may further expand the functions fin |. (77, u, e) as 



n^rfai^/w) = -/(^)[a n ,r(??,/")a r (/i,^) + /3 n ,r(^, A*)Sr(p, <p)], 



(4.13) 



where a n>r (r?,/i) and (3 air (r), fj.) are functions to be determined by the solutions to the 
Boltzmann equation, and 



dv 



(4.14) 



Now, plugging relation ( |4.13|) into ( |4.2j ), using the geodesic equation for the photon, 



dv 



dr\ 
d\' 



(4.15) 



where A is an affine parameter, and H is the Hubble parameter in conformal time, and 
expanding the photon momenta in spherical coordinates around the GW direction n, we 
obtain the Boltzmann equations for the radiative transfer in the presence of weak gravita- 
tional fields (dropping the indices n, r for the sake of simplicity) |3[ 

a 



33 



P(v, A*) + [q(v) + in-nMv, A*) 



16 



d_ 
dr) 



C(»7i/*) + [q(v) + «V*]f(»7, M) 



dr] 



Kv), 



where we have defined 



and 



£(77,//) = a(rj,n) +P(r),n), 



= J dp' 

and introduced scattering rate q(rj) defined by 

q(rj) = <j T N e (r))a{vi). 



(4.16) 
(4.17) 

(4.18) 

(4.19) 

(4.20) 



The solutions to the Boltzmann equations fl4.16|) and ( 4.17] ), given by the functions 
a(rj,/j,) and /3(?y, /x), are the essential elements for the computation of the anisotropies and 
polarization of the CMB, as we sketch in the next section. 

5. Harmonic analysis on a 2-sphere 

To begin with, let us first construct the polarization tensor associated with the Stokes 
parameters Q(8,(p) and U(6,ip), where the coordinates (#,</?) describe the position of a 
given region of the sky. The Stokes parameters Q and U can be cast into the symmetric 



trace- free (STF) polarization tensor 31f| , |44j] 



ab( 



Q -Usin6 
-Usin9 -Q sin 2 9 



(5.1) 



On the two-sphere tensor analysis can be easily implemented; the "divergence" and 
"curl" of a symmetric rank-2 tensors are respectively given by T ab , ab and T ab ]ac e c b , where 
":" denotes covariant differentiation, g ab and e ab are respectively the two-dimensional metric 
and antisymmetric tensors on the 2-sphere, given by 




£ab(0,¥) 



sin 




(5.2) 



With these elements on hand, we introduce invariants which can be built up from the 
polarization tensor V a b and its derivatives. From the symmetric tensor on the 2-sphere we 
construct two of the invariants, and from the second derivatives we construct the other two 
in the form of a "divergence" and a "curl" [33], 



9 ^ab, 



V 



ie ab P ab , 



E 



-2V, 



ab 



:ab 



B 



-2V ab - bc e a c . 



(5.3) 



EC 

o 
o 

o 
o 
-J 

o 
o 
o 



-14- 



With these invariants we get a very convenient way to completely characterize the radiation 
beam, since they do not depend on the reference frame chosen. We now proceed to expand 
the invariants (I,E,B,V) in spherical harmonics in order to perform an analysis on the 
each multipole of the radiation field p3| 



1(8, p) = Yl aJ m Y em (e,ip), 



=0 m= 



£=2 m=- 
oo I 



=2 m= 




X£-2)\ 



B 
(in 



Yi m (0,<p), 



(5.4) 
(5.5) 
(5.6) 
(5.7) 



=0 m= 



It is important to stress that these expansions are consistent with the similar definitions 



in the literature [31|, 32]. 

We are now in position to write down the I, E and B functions (|5.4Q , (|5.5| ) and (|5.€|) 
in terms of the functions a(fi,(p) and /3(fi,(p) introduced in ( |4.13|) . From (f4.1|) we obtain 
for a monochromatic radiation beam 



I(rj, v, 0,(p) = — [m(rj, u, 8, ip) + n 2 (v, v, 6, ip)} , 
Q(V, 8, ¥>) = — [(™i(?7, 8, P) ~ ^(rj, v, <p)\ 



U(ri,v,0 t <p) 



hv ; 



(5.8) 



so that from (4.12), Q4.13 ) and ( |5.8| ) we get (restoring the n-dependence of the Fourier 
expansion) , 

hu 3 



I nir (r,,u,9,ip) = ^[n(°V) + /(^„, r (w)(l-^)^]. 
Q^v^tp) = ^/(^)/3 n , r (r ? , / x)(l + / i 2 ) e ±2 ^, 



±2iip 



(5.9) 
(5.10) 
(5.11) 



From equations ( |5.9| - [5.1l| ) we may readily evaluate the expressions for /, E and B, 
using ( |5.1| ), (|5.2|), ( ^ ) and ( p.9[]5.11| ); then, integrating over photon frequencies, we obtain 



In,r (M) <P) 
E n ,r (fJ>, V) 



7 [(l- / u 2 ) «n,, (r?,//) e ±2 ^] 



"7 



<i 2 d 

l-/i 2 ) ( (1 + /i 2 ) _ + 8/i— + 12 ) /? n , r (r/,/i) 



Bn,r (H, <fi) = =F7 



2 (1 - M 2 ) ( + 4i|- ) /V fa, M ) e ±2i ^ 



(5.12) 
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where we have defined 

7=/ dv^rf{v). (5.13) 

Once we have the key expressions for evaluating the power spectrum correlation func- 
tion all we must do now is solving the Boltzmann equations ( |4.16 ) and ( f4,17| ), which we 
handle in the next section. 

6. The solutions to the Boltzmann equations 



In the paper |33[] the authors discuss an analytical method for solving the Volterra equation 
represented by ( 4,16| ) in terms of a series expansion, and compare their results with the 



exact numerical solutions. Here we follow only their numerical approach, which we sketch 
below. To do so, we introduce first the functions 

*M = j^s(i)i(v), (6.i) 

m = e -*W>m, ( 6 . 2) 

where the function r(ry) represents the optical depth of the universe, and is defined within 
a time interval rj' and r/: 

t( V , V ') = I dr/'W); (6.3) 
Jrf 

g(rj) is the visibility function, written as 

g {n) = q(T])e~ T(r,) = ^e- T ^. (6.4) 

Taking rj = 770, we further write the optical depth from a given conformal instant 77 to 
the present as T{rjQ,rj) = t(t]), that is 

rvo 

T ( 7 l) = / drj'q(rj'). (6-5) 

Jrj 



Now, using these definitions, the formal solutions to the equations ( 4.16j ) and ( 4.17| ) 



are given by the integral relations [33] 



p{n,n) = e r[ ^- inm r drf $(r/)e iW , (6.6) 







£(77, n) = e rM - intMri f V dr,' H{r,')e in ^' . (6.7) 







Now, since the function H(r/) is known, we can obtain a single integral equation for 
the function (|67T|) by plugging (|6^ ) and ( |6J|) into (|4.19| ), so that 



Zfa) = e^> jT dfulrf^ (1 + ^ 2 ) 2 <%') - i (1 - m 2 ) 2 F^)}^'^; (6.8) 



such expression can be further simplified by introducing the kernels K±(j] — 77'), defined as 



K±(v~rf) 
yields 



dfi(l ± ji A ) z e 



2\2 in/x(jj-j/) 



(6.9) 



(6.10) 



so that expression 

An) = e T{v) £ drfU+ir} - V 'Mv') - - r/)H(r/)\. 

The final equation for <£(??) is obtained by multiplying both sides of this equality by 
(3/16)q(rj)e~ T ^ and using the expression (6.1), so that 

Hv) = ^q(v) jf dr/^(rf)K+(n - rf) + G( V ), (6.11) 

where G(rj) is related to the function fl6.2| ), 

G{v) = ~q(fj) £ dr/H(rf)K.(r,-rf). (6.12) 

The solution to Volterra integral equation ( |6 . 1 1[) provides the values of the functions 
a and /3 for every conformal instant 77; in particular, to the present-day 770, the expressions 
a(r)o,fj,) = a(fi) and f3(r]o,fj,) = /?(//) are respectively given by 



no 



drj (H n>r (ri) - ® n ,r(v)) e 



-%nC 



no 



drj <& n ,r{rj)e 



-ifiC 



(6.13) 



where we have introduced the variable £ = n(rjo — 77). From ( |6.13] ) we compute the co- 
efficients af m : to do that we substitute expressions ( |5.4f|5.6| ) and ( |6.13 ) into ( |5.12| ), and 
integrate over angular variables, so that 



T 
n E 

u 'lm,nr 



i) 2 (<^2,m^l,r + 5-2,m52,r) af } nr' 
2 fx. x , x . x. \ „E 



-2,m02,s J a e , 



B 



where 



E 

a Lnr = 7 



JJ 

a £,nr 



rvo 

7x/4vr(2£ + l) / d V (H n , r ( V ) - $ n , r (v)) 2/(0, 
J 

rvo 

V^(2£ + 1) / drj $^(77)^(0, 
J 

rvo 

7VM2I + 1) / cfy $ nir (77)^(C), 
JO 



(6.14) 

(6.15) 
(6.16) 
(6.17) 



and ?XC)j ^(C)> ^t(C) are the multipole projection functions which appear after the 
integration over the angular variables, whose form are given by 



2KC) 

2^(0 

5/(0 



(£+2)1 MO 

(£-2)! C 2 ' 

2-^^)^(0-^-1(0 



(£-1) 



(6.18) 



tool of CMB physics: the correlation function Cf x ', where X, X' = E,B. Since each af n 
in ( |5,5| ), (p^), depends upon the wavenumber n and polarization state r, we write af m as 



Once we have obtained the expressions for |uj ) we can evaluate the most important 

x 

f.rn 

tm 

af m nr , so that the correlation function is given by : 
C 2 f e 

~ ^2(21 -\- 1) J n< ^ n ^ ] ^ ] [ fl fm,nr a fm,nr + a fm,nr a to,nr] ■ (6.19) 

T — 1 ; 2 711, = — ^ 

To conclude this section, let us say some words about the ionization history of the 
universe. We focus here on the time of decoupling, were the CMB radiation was released 
from the primordial nuclei. At decoupling, which took place around redshift z ~ 1088, the 
universe underwent a transition from a completely ionized state to a state in which neutral 
hydrogen and helium atoms were formed. In this process the radiation decoupled from 
the matter, originating the CMB radiation and a neutral pre-galactic baryonic medium 
(PGM). Then, at some redshift between 14 < z < 6 the PGM was ionized again by the UV 
radiation from the first luminous objects, leaving the intergalactic medium (IGM) ionized 
|45[| . Such process is called reionization, and would leave observable imprints on the CMB 
polarization spectrum due to the interactions of the CMB photons with the free electrons 
now available due to the reionized medium ^(|, [[I?]]. However, the reionization epoch is 
still not fully understood, and many models have been proposed to shed a light on the 
physics of this process (see j|8| and references therein), which can be homogeneous models 
with a sudden reionization (e. g. as discussed in |49|], f50|), or extended models with double 
reionization [51|, among others (see |4{^] for a more comprehensive list of papers). 



In the present paper we shall consider solely the epoch of recombination, whose physical 
process is very well understood. Despite reionization is fundamental to understand the low- 
multipole behavior of CMB polarization, it can be neglected in a first-approximation to 
study temperature anisotropies generated by the tensor modes. Next, as for the decoupling 
epoch, the formulae for the density of free electrons N e (rj) and the fraction of ionized 
electrons are given by [M] 



iV (77) = ( 1 — — ^ Xe ^ nbPc ( a ( 7? °) 



3 



and 



2 J m p \ a(rj) 



Yr,\ / c? \ / m„ \ „„„_i / z \c2-l / a' 



X J v )=[l-tP\ (—?—) — Of 1-1 ( — ) " [-)(l + z)-\ (6.20) 

yu v 2 J ^1000/ \2o t Rhpc) b V1000/ V«y 

In these formulae Y p ~ 0.23 is the primordial helium mass fraction, Qf, is the baryon content, 
and trip is the mass of a proton. The constants are given by c\ = 0.43, C2 = 16 + 1.8 In f2#; 
we take fl b = 0.046 |4]. 

7. CMB Anisotropies in Massive Gravity 



Let us now compute the Boltzmann equations for radiative transfer in the presence of weak 
gravitational fields in Massive Gravity. The anisotropy and polarization power spectra, as 



discussed in sections ||, || and |(| is a combination of two physical processes appearing in 
the Boltzmann equation ( |4,2| ): Thomson scattering, represented by the collisional integral 
fl4.3| ) on its right-hand side, and the gravitational redshift of the photon, represented by the 
geodesic equation ( |4.15| ) on its left-hand side. The collisional integral for Massive Gravity 
is the same as in GR, since the process of scattering involves only photons and electrons, 
and the lagrangian of matter is minimally coupled to the metric. The basic change is due to 
the different gravitational field strength represented by hij in the geodesic equation, which 
now depends on the details of the underlying gravitational model. However, the general 
form of Boltzmann equations ( |4.16| ) and ( 4.17|) are the same for both GR and Massive 
Gravity (see |2^] for details), which we proceed to evaluate numerically in what follows. 

We start our analysis of the possible signatures of massive gravitons on CMB anisotropies 
in the light of our discussion in section 3.2 and the results reviewed in section ||. For the 
sake of simplicity we consider first the case of GR, and then extend our discussion to Mas- 
sive Gravity. From figure || we see that long wavelength tensor modes (that is, low n, in 
the figure corresponding to n = 10) remain frozen during decoupling, so that they won't 
contribute significantly to CMB anisotropies and polarization at this epoch; this effect will 
be mainly due to the short wavelength tensor modes (high n, in the figure corresponding 
to n = 10 2 and n = 10 3 ; notice that tensor modes with n = 10 4 are "dead" at decoupling), 
which fell into the horizon at an earlier time. For late times the situation is inverse: long 
wavelength modes are now dead, and short wavelength modes are falling into the horizon. 
Therefore, this quick and simple analysis leads us to conclude that short wavelength tensor 
modes will be predominant over early times, and long wavelength modes at late times. 
In terms of multipoles of the radiation this feature can be understood in the following 
way: short wavelength modes and early times correspond to larger values of the variable 
C = n(r]Q — 7]), so that the spherical Bessel functions appearing in ( 6,18| ) will have nonzero 
values only for higher multipoles; conversely, long wavelengths and late times will give rise 
to smaller £, and hence je(C) will be nonzero only for low multipoles. 

Therefore, roughly speaking, in GR long wavelength tensor modes will influence mainly 
the CMB low multipoles, whereas short wavelength tensor modes will substantially con- 
tribute to the polarization on higher multipoles. This is also true in the case of Massive 
Gravity, but now with one substantial difference: long wavelength modes fall into horizon 
earlier than their massless counterparts, as seen in figures |2| and [|, altering then the form 
of the signature at low multipoles. Also, the oscillations noticed in the late-time evolution 
of the massive modes, valid not only for those possessing long wavelengths, but for the ones 
possessing short wavelengths as well (as seen in figure ||, for instance) would contribute to 
yield a distinct signature on low multipoles. 

However, since we are not considering reionization in this paper, the distinct signatures 
discussed above will not show up on the polarization spectrum. This happens because in 
this case the visibility function is zero at this epoch, and then the source functions ( p.l| ) 
will be zero, since they depend linearly on g(rj). However, this fact does not affect CMB 
anisotropies, since the mode coefficients aj m , given by ( |6 .15 ), depends on H(r]), defined by 
( |6.2| ), which is not zero even in the absence of reionization. Since the source function H(r/) 
depends on the tensor mode amplitudes, and massive and massless modes are different at 



late times, we conclude that the massive modes would leave a distinct signature on CMB 
low multipoles even without reionization. 

Let us now perform a numerical analysis in order to check the points discussed above. 
As we have discussed in Section 3.2, if the mass of the tensor mode is less or equal than 
mi, Massive Gravity produces the same results as GR; we choose then masses within the 
range m = 10~ 27 cm~ 1 - m = 10~ 26 cm -1 , whose associated anisotropies power spectrum is 
depicted in figure ||. Figure 1C show the low-multipole region of the correlation function ^ 
in detail. 



CMB Anisotropy for tensor modes 




I 



Figure 9: The correlation functions Cj for GR and Massive Gravity. Notice that the massive 
gravitons leave a signature on the spectrum for low multipoles. 



This figure show distinct signatures for massless and massive gravitons, as we have 
argued above. Therefore, for the range of masses selected, massive tensor modes leave a 
clear signature on low multipoles I < 30. Since the heavier modes fall into the horizon 
earlier, they have the stronger signature, as shown. If we had chosen a different mass, say 
m = W cm , the signature would be stronger, and possibly would appear for multi- 
poles t > 30. This can be explained by simply analyzing the trend shown in figures || 
and ||: the heavier the mass, the earlier the modes fall into the horizon, which correspond 
to higher multipoles. However, even in this the trend shown in figure [l(] indi- 

cates, the signature will be particularly strong on low multipoles. Therefore, if the tensor 
modes of the metric fluctuations are massive, they could be detected directly by the CMB 
anisotropy power spectrum if their mass are greater than the limit mi ~ 10~ 29 cm -1 , and 
their signatures would be noticeable specially on low multipoles. 

Therefore, the results above indicate clearly that the future measurements on the TT 
correlation might be decisive for probing the existence of massive tensor modes, for the 
signature left by them could be strong enough to be distinguished from those of the massless 
modes. 
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Figure 10: The low-multipolc "tail" in the TT correlation function. Notice the quite distinct 
signatures for i < 30 for the mass range selected. 

8. Conclusions 

In this work we have first studied the time evolution of massive tensor modes and shown 
that there is a graviton mass limit, mi ~ 1CP 29 cm -1 , such that gravitons with masses 
m <m>i behave indistinguishably from massless gravitons. The same happens to gravitons 
with short wavelengths (wavenumbers n > 100 in our example): their behavior is almost 
the same as of the massless gravitons for all the masses taken into account here. 

We have also shown that long wavelength massive tensor modes fall into the horizon 
earlier their massless counterparts, whereas short wavelength modes behaves quite similarly 
as in GR. The net effect of this behavior, as we have shown in the TT correlation function 
plotted in figures | and 0, is a distinguished signature on low multipoles; the heavier the 
mass of the mode, the stronger is its signature compared to that of massless gravitons. For 
the range of masses considered here, m = 10 -27 cm" 1 - m = 10" 2& cm" 1 , the signatures 
show up at £ < 30; however, we have argued that such signatures might appear at I > 30 
in the case of masses greater than m = 10 _25 cm _1 . 

Therefore, our results indicate that the future precise measurements of the CMB 
anisotropies induced by tensor modes might be decisive for probing the existence of massive 
gravitons, for the signature left by them could be strong enough to be distinguished from 
those of the massless modes. 
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